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Abstract 

Three-body recombination rates for cold 4 He are calculated with a new method which exploits 
the simple relationship between the imaginary part of the atom-dimer elastic scattering phase shift 
and the S-matrix for recombination. The elastic phase shifts are computed above breakup threshold 
by solving a three-body Faddeev equation in momentum space with inputs based on a variety of 
modern atom-atom potentials. Recombination coefficients for the HFD-B3-FCII potential agree 
very well with the only previously published results. Since the elastic scattering and recombination 
processes for 4 He are governed by "Efimov physics" , they depend on universal functions of a scaling 
variable. The newly computed recombination coefficients for potentials other than HFD-B3-FCII 
make it possible to determine these universal functions for the first time. 

PACS numbers: 21.45. +v,34.50.-s,03.75.Nt 
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I. INTRODUCTION 



Much of the current interest in the properties of cold dilute atomic systems arises from the 
possibility of examining exotic three-body effects which are generically referred to as "Efimov 
physics" . Such effects arise when the two-atom (dimer) scattering length a is large compared 
to I the range of the atom-atom potential which is typically of order the atomic size. Some 
remarkable "universal" features emerge in this limit (see, e.g., refs. p], 0] and references 
therein). For example, it is rigorously true that, as the dimer scattering length diverges 
and the binding energy goes to zero, the atom-dimer scattering length also diverges and the 
three-atom (trimer) system develops an infinite tower of bound states whose energies form 
a geometric sequence with zero as an accumulation point. In physical systems, departures 
from this limiting case scenario typically arise when, for example, If a is finite. However, 
when such departures are small, universal Efimov physics will still dominate. 

The history of Efimov physics resembles to some extent that of Bose-Eistein Condensation 
(BEC) in that the theory was well in hand long before experimenters were able to confirm its 
predictions. Indeed, it was only within the past year that the Innsbruck group , exploiting 
the technique of Feshbach resonances in an ultracold gas of 133 Cs atoms, saw evidence of 
Efimov physics in three-body recombination rates. Most of the elements of their analysis of 
the data are summarized in a recent review article by Braaten and Hammer 0]. Such an 
analysis, based on concepts of universality, requires some independent empirical input. A 
recent paper by Braaten, Kang and Platter [4| examines the universality constraints available 
from cold 4 He and how these apply to the case of 133 Cs. 4 He is an especially interesting 
case because (i) it is relatively simple theoretically and (ii) nature has kindly provided a 
dimer scattering length which is large enough (~10 nm) that the criterion for Efimov physics 
to appear is generously satisfied. On the theoretical side, atom-atom potentials have been 
developed to a high degree of sophistication. Four widely used potentials are HFD-B 
HFD-B3-FCII @, LM2M2 Q and TTY [8]. These have been employed in numerous "ab 
initio" calculations of dimer and trimer properties (see, e.g., refs. [9|, llO, llll ) • These results 
in turn have been used as input for the universality-based analyses mentioned above. On 



the experimental side, 4 He dimers have been observed 12] , Il3| and the atom-atom bond 



length has been determined to be 5.2±0.4 nm from which an atom-atom scattering length 
of 10.4l°;g nm has been deduced. (Scattering lengths obtained using the four atom-atom 
potentials mentioned above are consistent with this range of values.) 4 He trimers in their 



ground state have also been observed [13|, [14J and a recent paper [15[ reports a bond length 



of 1.1 

-a5 nm which is consistent with theoretical estimates. This same paper indicates that 
the trimer excited state was not observed. The upper limit for the concentration of excited 
trimer states was determined to be appreciably less than theoretical estimates, a result which 
led the authors to entertain the possibility that such a state does not exist. 

It is clear from the above discussion that three-body recombination rates are crucial in 
establishing the extent to which Efimov physics can be realized in the laboratory. To date, 
however, only a single "ab initio" calculation exists. It was performed by Suno et al [l6[ for 
4 He using HFD-B3-FCII. As pointed out in Ref. 0], this is unfortunate because it greatly 
inhibits the unversality-based analysis. In the present work, we present new results for 
the 4 He three-body recombination coefficients. These quantities are obtained for the HFD- 
B, LM2M2 and TTY potentials as well as for HFD-B3-FCII where excellent agreement 
with Suno et al. is achieved. As will be discussed below, a new calculational technique is 
employed which appears to be accurate and computationally efficient. Finally, we employ 
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HFD-B3-FCII 


1.58730 
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2.624 [17J 


0.2835 


12.017 


7.06 xlO" 28 


8.738 


1.1477 


HFD-B 


1.68541 


8.856 


2.74 [10] 


0.2872 


12.149 


4.10 xl0~ 28 


8.480 


1.1222 


LM2M2 


1.30348 


10.018 


2.28 [10] 


0.2932 


11.422 


24.7 xl0~ 28 


9.643 


1.2367 


TTY 


1.30962 


9.996 


2.28 [10] 


0.2855 


11.548 


22.7 xl0~ 28 


9.620 


1.2270 



TABLE I: 4 He dimer and trimer properties. Trimer results are from FF calculations. 
= K 3 for i?fe rea fc M p — > 0. See text for definition of quantities shown. 



our new recombination coefficients to extract for the first time the universal functions, hi 
and fi2, introduced in Ref. [H which govern low breakup energy recombination rates when 
Efimov physics is dominant. 



II. CALCULATION OF RECOMBINATION RATES 

The rate equation for the density of a cold thermal atomic gas is given by [1 



dn A 3K 3 3 

— 7— = n-A 1 

dt 6 A K ' 

where n A is the atom number density, K 3 is the recombination coefficient and where we have 
assumed that all three atoms involved in the recombination leave the trap. We also have 
ignored the contribution from collision induced dissociation processes. For 4 He atom-dimer 
scattering, the only inelastic channel corresponds to three-body breakup. Therefore the 
recombination ^-matrix is directly related to the inelastic part of the atom-dimer S-matrix. 
It follows 3 that 



K 3 = ^(2L + 1) K 3 (L) = £(2L + l^lf (2) 

L L 

where fi = rriHe/ and k depends on the energy, E, in the breakup channel via k = \f2jlE. 
Also, L is the angular momentum in the scattering channel and S™ is the inelastic S-matrix 
element in that partial wave. Via the argument presented above, we have 

\S™ | 2 = 1 -e" 4 ^ (L) (3) 

where Si(L) is the imaginary part of the phase shift for elastic atom-dimer scattering. Hence, 
calculation of K 3 (L) is equivalent to finding the atom-dimer elastic phase shift above the 
breakup threshold. We accomplish this by solving the momentum-space Faddeev equation 
for three-body scattering, to which we now turn. 

Following, e.g., Watson and Nuttall [HI], we assume the s-wave atom-atom interaction is 
described by a separable potential: 

V(q,q')=g(q) X g*(q') (4) 

where g(q) is a form factor (FF) to be discussed at length below. The atom-atom t-matrix 
is then 

t{q,q';E + )=g{q)r{E + )g{q') (5) 



3 




FIG. 1: (Color online). Atom-dimer phase shifts using potential LM2M2. Circles are from Mo 



tovilov et al. [lOj]; crosses are from Roudnev 11]. To be compared with Figure 3 of Ref. 18]. 
Breakup threshold is at E cm + B 2 = 1.30348 mK. ET-ERE NLO and NNLO calculations are 
virtually identical to the FF curve on this plot. 



where E + = E + ir], r(E + ) = A/[l - A h(E + )} and 

h(F+] f d 3 k \g(k)f 

where m is the atomic mass. Following, e.g., the standard Faddeev-based treatment of 
Ref. [Hj], we nn d that the fully off-shell atom-dimer t-matrix for three identical bosons in 
the L-th partial wave is 

T L { P ,k-E + ) = 2Z L { P ,k-E + ) + 2 J -0- 3 Z L (p,q;E + )r^E + -^j T L (q,k;E+) (7) 

with 

1 f +1 i n r g(p/2 + k) g'(p + k/2) 



where Pl{x) is a Legendre polynomial. (The full derivation of Eq. [7] is somewhat lengthy. 
Ever longer is the full history leading to the development presented in Ref. 19|. Neither the 
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derivation nor the history can be fully recounted here but suffice it to say that the many 
references presented in Watson and Nuttall pi| dating from the 1950's and 1960's contain 



the direct antecedents to, e.g., the trimer bound state calculations of Lim et al. [20|. The 



expressions for the three-body amplitudes contained in these early references have been re- 



derived many times since; Refs. [2l|, |22|, |23|, |24], |25| are but a few representative examples. 



Note that, in virtually all of the these more modern papers, two-body contact interactions are 
assumed. See, however, Afnan and Phillips [26[ for a modern reference in which expressions 
very similar to ours appear even though, in their work, the form factors g(q) are introduced 
for pedagogical purposes only and are set to 1 early on.) The elastic scattering amplitude 
is then given by 

f(0) = Y,(2L + l) f L (p) P L (cos6) (9) 

L 

where 

HP) = ~ T L (p,p; ^-B 2 + ir))= (10) 
Sti 4m Zip 

where the atom-dimer t- matrix is fully on-shell. Here B 2 is the dimer binding energy which 
satisfies 1 — Xh(—B 2 ) = and the bound state normalization Z is given by 



z -i/2 = dh ( E ) 



olE 



(11) 



-B 2 



If the form factor g(q) is taken to be a 5-function, then h(E + ) = rn ^-^ lE+ after the di- 
vergence is removed. Setting Att/Xtu = a 2 , the atom-atom scattering length, we find that 
Eq. [7J is the standard Effective Theory (ET) result at leading order (LO) in the Effective 



Range Expansion (ERE) of the atom-atom t-matrix. (See, e.g., Refs. [18l. l27j for extensive 
discussions of this treatment and for additional related references. We will subsequently 
refer to results based on this formulation as ET-ERE calculations.) Below the breakup 
threshold, the resulting Faddeev equation, Eq. [7J can readily be solved as it stands by dis- 
cretizing in momentum space. This yields a system of simultaneous linear equations which 
can be handled with standard linear algebra packages such as LAPACK. Above breakup, 
as noted long ago by Aaron and Amado [28[ in the context of neutron-deuteron scattering, 
Zl becomes ill-behaved along the real axis. A simple and efficient solution to this prob- 
lem was given by Hetherington and Schick [29| even longer ago. They suggested deforming 
the contour of integration over q in Eq. [7J away from the real axis, solving the Faddeev 
equation for Ti of complex argument, and then using the Faddeev equation once more to 
obtain Tl (p, p; — B 2 + irj) . Finally, the relation given in Eq. [10] is utilized to extract 
the phase shift, Sl(p). When carefully applied, this prescription can yield accurate complex 
phase shifts above the breakup threshold and hence, via Eqs. [2] and [3l numerical values 
for the recombination coefficient, K 3 . It should be noted that care is required because the 
imaginary parts of the phase shifts are typically quite small. For L = near threshold, Im 
5l is typically of order 10~ 5 degree! The fact that we compute values of Im b~L which are 
typically less than 10~ 8 degree just below threshold is an important necessary condition for 
trusting the small values we find above threshold. Another indication that our numerics are 
adequate in computing Im 5l is that our values of K 3 (L = 0) go to a constant as threshold 
is approached, as required by general considerations [30|, 3l[. This requires considerable 



numerical precision since, in this region, Im 5q is proportional to E 2 , where E the total 
energy in the breakup channel. 
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Ref. [10] 


Ref. [32] 


Ref. [11] 


Ref. [33] 


Ref. [34] 


Ref. [18] 


Present work (FF) 


HFD-B 


13.5±0.5 




12.19±0.01 




12.47 




12.149 


LM2M2 


13.1±0.5 


12.6 


11.54±0.01 


11.425 


11.87 




11.422 


TTY 






11.58±0.01 




11.89 


11.573 


11.548 



TABLE II: 4 He atom-dimer scattering lengths in nm compared with previous results. 



In the following section we present results using Eqs. [7]and[T0l In addition, we have done 
calculations at next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) in 
the ET-ERE 18|, |27j. To facilitate comparison with the published results of the latter ref- 
erence, we present our definitions of these calculations. Recall the effective range expansion 
of the two-body s-wave scattering amplitude: 



p cot S 



1 

a 2 



-rp 2 + 0\p% 



(12) 



where a 2 is the scattering length and r is the range parameter. We note that for 4 He atom- 



■V 



atom scattering, a 2 — 10 nm while r ~ 0.75 nm. Then, at LO, the quantity r(E^ 
Eq. [7] may be expressed as 



m 



T 



4m 



An 



1 
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-7 + V / !(? 2 -^)+7 2 



(13) 



where 7 = Atc/ Am and where we have used 
At higher orders in the ET-ERE, we have 



k 2 — / y + irj)/m. At LO, we take 7 = l/a 2 . 



r [E + - 
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47T 
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where N n 



1 or 2 for NLO or NNLO, respectively. In this expression, we take 



(14) 
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in 



2r 
1 

a 2 



(15) 



which means that p = 27 gives the position of the bound state pole in the two-body T-matrix 
which corresponds to a zero in 



V/o(p) =pcotS -ip 



(16) 



when terms of 0\p 4 } are dropped in the effective range expansion for pcot So (Eq. [T2|) . Note 
that the expression for r(E + — in Eq. [H] ensures that the pole structure of that quantity 
is the same at NLO and NNLOls at LO [23]. 

In addition to ET-ERE calculations, we will also show results for which the form factor, 
i.e., the g(q) appearing in Eqs. H] and El is not assumed to be a delta function but rather is 
taken to be of the form 



9(o) 



1 + 6 P\ 



[1 + g 2 //? 2 ) 2 



(17) 
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where {3 and C are fixed by the low energy atom-atom phase shifts [35(. We compute these 
phase shifts directly from the atom-atom potentials along with the corresponding dimer 
binding energies. 

The dimer binding energies and atom-atom scattering lengths for the four potentials 



considered here appear in Tabled! As discussed in detail in, e.g., Refs. [la. l2ll. |22|. |23|. 136 
the Faddeev equation for L = is a special case of considerable interest as it is here that 
Efimov physics is manifest. In the standard ET-ERE approach, T^ =0 at LO exhibits a log- 
periodic dependence on the upper limit, A, of the integration over q in Eq. [7J Moreover, a 
three-body force enters at lowest order in the calculation and is accounted for via 

Z L=0 -> Z L=0 + ^ (18) 

where the three-body interaction strength Hq is dimensionless. For the ET-ERE calcula- 
tions reported here, we choose A arbitrarily apart from the requirement that it to be large 
enough that the overall results do not depend on its magnitude. Then H is adjusted to 
reproduce B% , the binding energy of the trimer excited state which is known to be of an 
Efimov character [2(|. (Note that our trimer binding energies are defined so that zero corre- 
sponds to the three particles being at rest and far enough apart that they do not interact.) 
Consequently, the value of H is A-dependent. When using the finite-range form-factors, the 
situation is somewhat different. Because the form-factors regulate the integrand in Eq. [7J 
there is no log-periodic dependence on A. The three-body interaction must still be included, 
with its value determined as before. Now, however, H is independent of the momentum 
cutoff. A thorough discussion of form-factor calculations like these will be presented in a 
future publication j35|. Values of B% and Hq for form- factor calculations (FF) are given in 
Table [J for the four atom-atom potentials we consider. 



III. RESULTS 



To begin this section, we focus on results for the atom-dimer scattering length. Table 
HT1 compares our FF values with those of previous works. Agreement with Roudnev 11] , 
Pen'kov 33| and Platter and Phillips [18] is quite good. The former calculations are based 



on differential Faddeev equations in configuration space in a total angular momentum repre- 
sentation. The author claims a two order of magnitude improvement in accuracy over earlier 
"ab initio" calculations of a similar nature by Motovilov et aL_ 13 
results by the latter group appear in Table [Til under Ref. 



34 



Note that more recent 
Platter and Phillips 
report a version of ET-ERE calculations using TTY. We find especially close agreement 



for LM2M2 with Pen'kov [33J] whose work we became aware of after the present work was 
completed and whose method is similar in some respects to ours. We are able to confirm 
an observation by Pen'kov that p cot So for 4 He atom-dimer scattering has an unusual low- 
energy behavior which greatly restricts the range over which the low order effective range 
expansion is accurate. We also note that this behavior necessitates some care in finding the 
scattering length by, in effect, extrapolating pcot<5 to p = 0. In particular, extrapolations 
which do not take into account this unusual behavior will yield scattering lengths which 
are too large. Such problems may account in part for some of the larger discrepancies with 
previous results appearing in Table [UJ Figure [^compares our atom-dimer L = phase 



shifts for LM2M2 with those of Motovilov et al. [lOj and of Roudnev [llj . We show results 



of our ET-ERE calculations at LO, NLO and NNLO along with form factor calculations 
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FIG. 2: (Color online). K^(L) computed using the form factor g(q) (FF) and atom-atom potential 
HFD-B3-FCII. Present results (dashed) are compared with those of Suno et al. 16f] (solid). 



(FF). (We remind the reader to look at Refs. [la , l27j ] for thorough discussions of ET-ERE 
calculations.) First we observe that there are significant differences between LO and higher 
order ET-ERE results but the ET-ERE calculations are well-converged at NLO. Moreover, 
the NLO and NNLO ET-ERE results agree well with the FF phase shifts. As will be dis- 
cussed in a future publication (35|, this behavior is readily explained. First, range effects are 
small here since r/a 2 for 4 He atom-atom scattering (a 2 and r are the two atom scattering 
length, respectively, in the effective range expansion of pcot<5) is roughly 0.75 nm / 10 nm. 
Also, the FF calculations include, by construction, certain of the finite range effects to all 
orders (35|. Figure [2] shows the present phase shifts lie consistently below those of Roudnev 
but are in good agreement with the calculations of Motovilov et al. All our calculated phase 
shifts (except for ET-ERE LO) are indistinguishable from those of Platter and Phillips [18]. 



Our values for scattering lengths computed with different atom-atom potentials and using 
different methods of calculation are displayed in Table IIHI For the ET-ERE calculations 
at NLO and NNLO, our values depart from those of Platter and Phillips [l8| only at the 
level of +0.27 and -0.14 %, respectively. Based on the small value of r/a 2 — 0.075, we 
expect the ET-ERE calculations to converge rapidly, an expectation which is readily con- 
firmed by examination of Table HTT1 Moreover, it appears that the ET-ERE calculations are 
consistently converging toward the FF results. Such behavior would be consistent with the 
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Pot 


ET-ERE LO 


ET-ERE NLO 


ET-ERE NNLO 


FF 


HFD-B3-FCII 


9.833 (-18.2) 


12.029 (+0.10) 


12.024 (+0.06) 


12.017 


HFD-B 


9.957 (-18.0) 


12.163 (+0.12) 


12.160 (+0.09) 


12.149 


LM2M2 


9.269 (-18.8) 


11.446 (+0.21) 


11.435 (+0.12) 


11.422 


TTY 


9.388 (-18.7) 


11.566 (+0.16) 


11.557 (+0.08) 


11.548 



TABLE III: 4 He atom-dimer scattering lengths for various potentials and calculation types. Values 
in parentheses are percent deviations of ET-ERE calculations from the FF results. 




^breakup ^-M^] 

FIG. 3: (Color online). K^(L = 0)'s computed using the form factor g(q) (FF) for various atom- 
atom potentials. The solid lines represent full calculations while the dashed curves are computed 
using the universal functions h± and /12 and the values of 0,2 and ao* for potentials HFD-B3-FCII 
and TTY. Universal calculations for HFD-B and LM2M2 agree with their full counterparts by 
construction. See text for discussion. 



assertion made above that the FF calculations include some range effects to all orders [35 
Table ULTI also makes it clear that there are significant differences between the HFD-B3-FCII 
and HFD-B atom-atom potentials on the one hand and the LM2M2 and TTY potentials on 
the other. Such differences will reappear in the discussions to follow. 

We now turn to evaluation of the 4 He three-body recombination coefficients. FF results 
for K 3 (L = 0) in the low energy limit are displayed in Table [H Our value for HFD-B3-FCII 
of 7.06 x 10 -28 cm 6 /s is in striking agreement with the only previously published value of of 



9 



0.05 r — i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 
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FIG. 4: (Color online). Universal functions h\(x) and ti2(x) as denned in Ref. [4j] and appearing 
here in Eq. 1221 The circles (hi) and triangles (/12) are values extracted from HFD-B3-FCII and 
LM2M2 calculations of K 3 (L = 0). The fits are: h x (x) = 0.0179358 x -0.379918 x 2 + 0.22438 x 3 - 
0.161562 x 4 and h 2 {x) = 0.0214188 x - 0.152046 x 2 + 0.029073 x 3 + 0.051001 x 4 . 



this quantity, namely 7.09 x 10 cm /s as reported by Suno, Esry, Greene and Burke 16 



Figure [2] compares our results for K 3 (L), L < 3 as a function of breakup energy with those 
of Suno et ai. The agreement is again quite good (apart, perhaps, from the location of the 
minimum for L — 0) and we conclude that our method of calculating these quantities is 
trustworthy. (We note that the converged ET-ERE calculations of K 3 are nearly identical 
with the FF results in Table d and Figure [2j) Thus we can for the first time compare 
the HFD-B3-FCII for K3 results with those employing other atom-atom potentials. As is 
apparent in Table [J the HFD-B result for K%(L = 0) in the low energy limit is similar to 
the value for HFD-B3-FCTI. The values for LM2M2 and TTY are similar to one another 
but larger than the other two by roughly a factor of five. The same behavior is also evident 
in Figure [3] where FF results for K 3 (L = 0) are plotted as a function of breakup energy. 
Nielsen and Macek [30] as well as Esry, Greene and Burke [9[ found that the magnitude and 
shape of this quantity depends sensitively on an interference of two distinct amplitudes. It 
is this interference which, for example, gives rise to the so-called Stueckelberg minima at 
Ebreakup = 2000-3000 /iK. We speculate that, in the present work, the interference is between 
the purely two-body contributions and the three-body contribution appearing in Eq. [TH In 
any case, all calculations display Stueckelberg minima which are qualitatively similar to that 
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found originally in Ref. 30J and later in Ref. 16 



We next address extraction of universal functions which describe the recombination co- 
efficients in the scaling limit @, 3?J as outlined by Braaten, Kang, and Platter Q]. They 



find that, under a few physically reasonable simplifying assumptions, the recombination 
coefficient takes on the form (their Eq. 23) 

K 3 (L = 0) = C max \ sin[s ln(a 2 /a *)](l + h^x)) + cos[s hi(a 2 / a^)]h 2 (x)\ 2 hal/m (19) 

where the scaling variable is 

x = (ma 2 2 E breakup /h 2 ) 1/2 , (20) 

Cmax is a universal constant with a numerical value of about 402.7 and m is the atomic 
mass. The quantity a 2 depends on the dimer binding energy, B 2 , via 

a 2 = ^B 2 jh (21) 

while ao* is extracted from the value of K 3 (L = 0) in the E brea k up — > limit by solving Eq. 
10 of Ref. [3], namely 



- Q) - 76871 " 2 ( 47T ~ 3 v / 3) sin 2 [g ln(a 2 /a 0i ,)] fia\ 
sinh 2 (7rs ) + cos 2 [s ln(a 2 /a *)] m 



Values of a 2 and a 2 /a * for the four potentials we consider are given in Table [0 Our values 
for a * differ by less than 3% from those we obtain using the prescription of Ref. Q , namely 
ao* ~ 0.32 k' 1 where 



2/i li2 Bg ] /h (23) 

and where /ii. 2 is the atom-dimer reduced mass while B^J is the ET-ERE LO trimer excited 
state binding energy in the limit that a 2 —>■ +oo. We note that a value of a 2 /ao* — 1.15 for 
HFD-B3-FCII is reported in Ref. (4j which is in accord with our result of 1.1477. 

Because we have results for K%(L = 0) for several atom- atom potentials, we can extend 
the work in Ref. j3] by solving Eq. [19] for the two functions, hi(x) and h 2 (x), for the first 
time. We accomplish this by using the values of K^(L = 0) obtained with potentials LM2M2 
and HFD-B since these yield the highest and lowest values, respectively, of the L = 
recombination coefficient (see Table [Hand Figure [3]). Our determinations of hi and h 2 are 
plotted in Figure [H The range of the h^s extends up to x ~ 1.1 which is just below 
the Stuekelberg minima where the assumptions used to obtain Eq. [J9] are no longer well- 
satisfied 0| . We have made polynomial fits to hi and h 2 and the results are also shown in 
Figure HJ the explicit forms of the fitted polynomials appear in the caption of Figure [H We 
note that hi is much larger than h 2 and that both functions go to zero as x goes to zero as 
should be expected from Eqs. [22] and [T9l 

Figure [5] shows sums of i^3(L)'s for L < 3 in scaling units for a variety of values of a 2 
and a *. The results for a 2 /a * = 4.76, 1 and 1.1477 can be compared directly with Figure 
2 of Ref. [I| where only uncertainty bands appear. The present results are consistent with 
those of Ref. [IJ except for the case of a 2 /a^ = 1 where our values fall below their narrow 
band. For this value of a 2 /a *, only h 2 (x) contributes (see Eq. [L_H]) and, as mentioned above, 
we find h 2 (x) to be very small. It thus seems possible based on what appears in Figure [5] 
that, for systems where a 2 /ao* — 1, K 3 could be two or three orders of magnitude smaller 
than the already small values found for potential HFD-B over a range of Eb rea k up up to 
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FIG. 5: (Color online). Sum of Ks(Lys for L < 3 computed using the universal functions, h\{x) 
and h<2(x), shown in Figure E] and various vales of 0,2/0,0*- These results are to be compared with 
Figure 2 of Ref. 



100 /iK (or x < 0.2). We also show in Figure [5] calculations employing LM2M2 and TTY 
for which the values of 0,2 /ao* and K 3 are relatively large. Dimensionful counterparts of the 
HFD-B3-FCII and TTY curves calculated using the universal functions are compared with 
the corresponding full calculations in Figure Agreement is excellent and it is gratifying 
that the universal calculations can account for the dependence of K%(L = 0) over the full 
range of variability caused by differences among the four atom-atom potentials we consider. 
(Note that the agreement for HFD-B and LM2M2 is exact by construction.) 



IV. CONCLUSIONS 

We have performed calculations of three-body recombination rates for cold 4 He using 
a new method which exploits the simple relationship between the imaginary part of the 
atom-dimer elastic scattering phase shift and the S-matrix for the recombination process. 
We compute the elastic phase shifts by solving three-body Faddeev equations in momentum 
space. These equations are relatively simple under our assumption of separable atom-atom 
interactions 35(. When the form factors which characterize the separable interactions are 



taken to be 5-functions, our elastic scattering calculations are identical to Effective Theory 
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treatments based on the Effective Range Expansion (ET-ERE) at lowest order (LO) in the 
effective range expansion of the atom-atom potential. (See Ref. [27| for relevant discussion 
in the context of neutron-deuteron scattering and Ref. [18J where 4 He atom-dimer scattering 
is addressed.) We perform such calculations using four widely used atom- atom potentials, 
namely HFD-B3-FCII0, HFD-B [J, LM2M20 and TTyQ. We also report ET-ERE re- 
sults at next-to-leading order (NLO) and NNLO 27| in the expansion parameter rja^ where 
r is the range parameter in the effective range expansion of p cot Sq for atom-atom scattering. 
This quantity is quite small for 4 He (~ 0.075) and convergence is rapid. We also present 
results for a new kind of Faddeev calculation employing finite range form factors whose 
momentum dependence is fixed by fitting to pcot^o- These calculations have some similar- 
ities to neutron-deuteron scattering calculations done long ago, shortly after the Faddeev 
equations were formulated 19, 28[. Recent work by Pen'kov 33] is also similar to ours. The 



ET-ERE calculations for the atom-dimer scattering lengths appear to be converging toward 
the FF results. Our converged ET-ERE results as well as our FF calculations of atom-dimer 
scattering lengths and phase shifts agree well with rece ntly published values by Motovilov 



et al. 10], by Roudnev 111) and by Platter and Phillips [18]. 

We have calculated three-body recombination coefficients, K S (L), for L < 3 using the 
four atom-atom interactions we consider. The HFD-B3-FCII results agree very well with 
the only previously published values of Suno et al. [l6j]. Our calculations using the other 
three potentials reveal that the recombinations coefficients vary by as much as a factor of 
six depending on these inputs. More specifically, the results for HFD-B3-FCII and HFD-B 
lie in the lower end of the range while LM2M2 and TTY values are nearly identical and 
establish the upper end. 

To the extent that the dynamics of the 4 He atom-dimer system are governed by Efimov 
physics, it should be possible to describe them via scaling arguments 0, 37|. In a recent 



work [4], Braaten, Kang and Platter have shown that, under a few physically reasonable 
assumptions, the L = recombinations coefficient can be expressed in terms of "universal" 
quantities including two universal functions of a scaling variable. Because they only had 
access to the HFD-B3-FCII results of Suno et al. [16], they were unable to determine these 
functions with any precision. Having in hand recombination coefficients for four different 
potentials, we are able for the first time to determine these functions with considerable 
accuracy. Indeed, when they are fixed using the HFD-B and LM2M2 results, the universal 
relation describes the HFD-B3-FCII and TTY recombination coefficients very well. We 
also note that one of the universal functions we determine is quite small which suggests 
that extremely small recombination coefficients can result if certain physical quantities are 
appropriately tuned. 
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